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SUMMARY 

The  response  of  accelerometers,  with  second-order  filters,  is  considered 
2 

for  i sin  and  sin  acceleration  pulses.  The  time  history  of  their  integrated 
(velocity)  response  is  also  treated. 

Graphs  relating  the  error  in  measuring  peak  acceleration  to  the  filter 
parameters  are  given,  and  the  results  are  related  to  practical  situations 
experienced  when  calibrating  accelerometers  in  RAE. 

The  effect  of  accelerometer  non-linearity  when  calibrating  by  velocity- 
change  methods  is  briefly  investigated  and  a method  of  measuring  the  non- 
linearity is  suggested. 
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1 INTRODUCTION 

These  notes  are  intended  to  serve  as  an  aid  to  quantifying  some  of  the 

errors  which  arise  in  the  calibration  and  use  of  shock  accelerometers.  Calibra- 

1 2 

tion  errors  associated  with  an  RAE  airgun  facility  have  already  been  considered 
particularly  in  the  context  of  signal-conditioning  and  calibration  equipment. 
Nevertheless,  Ref  2 identified  important  sources  of  error  related  to  the 
characteristics  of  the  accelerometer  under  test  and  it  is  felt  that  some  of  thes>*» 
warrant  further  consideration. 

Accelerometers  for  shock  measurement  must,  in  general,  have  sufficient 
acceleration  range  and  frequency  response  to  register  the  shock  waveform 
adequately.  Thus,  most  accelerometers  used  for  measuring  shocks  of  up  to  say 
lO^g^  , are  of  the  piezoelectric  type  with  undamped  natural  frequencies  of  the 
order  of  100  kHz.  These,  being  undamped  , are  likely  to  ring  when  subjected  to 
shock  and  moreover,  other  resonances  associated  with  the  mounting  of  the  instru- 
ment are  likely  to  colour  the  output  waveform.  Such  resonances,  more  often  than 
not,  completely  obscure  the  shock  waveform  and  hence  it  is  usual  to  place  a 
filter  in  the  measurement  system  in  order  to  obtain  what  is  hoped  to  be  an 
adequately  faithful  representation  of  the  shock  pulse.  The  choice  of  filter  and 
where  to  put  it,  is  more  art  than  science. 

Filters  may  be  electrical  or  mechanical.  The  former  are  usually  employed 
since  they  are  easy  to  use  and  are  readily  adjustable  to  suit  the  prevailing 
conditions.  But  extreme  caution  must  be  exercised  lest  they  mask  overloading  of 
the  accelerometer,  or  its  conditioning  unit,  by  inputs  whose  spectrum  contains 
frequency  components  outside  the  pass  band  of  the  filter.  It  is  therefore  better, 
in  principle,  to  introduce  a mechanical  filter  at  the  very  beginning  of  the 
measurement  system.  Such  filters  are  available  commercially  in  the  form  of 
accelerometer  mounts,  but  their  use  is  limited  to  small-shock  measurements. 

Simple  calculations  quickly  reveal  the  problems  of  designing  such  filters  with 
useful  amounts  of  damping  to  withstand  shocks  of  10  gn-  Nevertheless,  it  is  the 
writer's  belief  that  such  filters  must  be  developed  and  used  if  shock  signatures 
are  to  be  more  reliably  measured. 

Because  the  characteristics  of  mechanical  filters  at  high  shock  levels  are 

unlikely  to  lend  themselves  to  accurate  control,  there  is  a need  to  quantify  the 

performance  of  accelerometers  with  such  filters.  For  reasons  of  simplicity,  the 

treatment  below  is  confined  to  second-order  filters  and  to  shock  pulses  of  4 sin 
2 

and  sin  form.  Many  shock  pulses,  and  particularly  those  generated  by  the 
airgun  used  for  shock  calibration,  are  roughly  of  these  forms. 


A special  case  of  a mechanical  filter  exists  when  calibrating  by  the  com- 
2 

parison  method  in  which  the  test  accelerometer  is  mounted  on  a reference  trans- 
ducer. Errors  can  arise  from  so-called  ’relative  motion'  between  the  sensing 

element  of  the  reference  transducer  and  the  surface  on  which  the  test  accelero- 

meter is  mounted.  The  coupling  between  the  reference  and  test  accelerometer  can 
be  represented  with  fair  accuracy  by  a second-order  mechanical  filter  whose  un- 
damped natural  frequency  and  damping  ratio  depend  on  the  design  of  the  reference 
transducer  and  the  mass  of  the  test  accelerometer.  This  case  is  analysed  in  the 
present  paper. 

For  completeness,  the  response  of  filtered  accelerometers  with  integrators 

2 

is  also  considered.  This  relates  to  calibration  by  the  velocity  method  in  which 
the  integrated  output  of  the  accelerometer  is  equated  to  the  measured  velocity 

change.  The  use  of  such  data  need  not,  in  fact,  be  confined  to  calibration  but 

may  be  extended  to  practical  measurement  situations  where  the  velocity  change  is 
known  or  can  be  guessed  with  fair  accuracy.  In  such  cases,  the  extent  to  which 
the  integrated  output  of  the  accelerometer  agrees  with  the  known  velocity  change 
can  provide  a valuable  indication  of  the  integrity  of  measurement. 

Whilst  emphasis  is  given  to  presenting  the  output-time  relationships  of 
accelerometers  and  deriving  the  associated  errors,  the  effects  of  accelerometer 
non-linearity  are  briefly  considered  in  the  context  of  velocity  calibration. 
Piezoelectric  accelerometers  tend  to  become  more  sensitive  at  high  levels  of 
acceleration  and  establishing  non-linearity  is  part  of  the  calibration  procedure. 


RESPONSE  TO 


SIN  PULSE 


Consider  a { sin  pulse  (Fig  la)  where  the  instantaneous  acceleration  is 
given  by 


a(t)  = a sin  iot 


0 < t £ If/u)  (1) 
t > it/u)  (la] 


where  a = peak  acceleration 


T = period  of  pulse 


time  . 


If  this  pulse  is  passed  through  a second-order  filter  of  transfer  function 
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where 

and 
it  is 


where 


s + 2c a)  s + u 
n n 

h(s)  * transfer  function  of  filter 

(o  *=  2-rrf 

n n 

fn  * undamped  natural  frequency 

C = damping  ratio  (assumed  < 1) 

shown  in  Appendix  A that  the  output  of  the  accelerometer  and  filter  is 


y(t)  * g(t)  + g(t  - J)H(t  - J) 
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whence 


lim  v(t) 

t-H» 


lim  sv(s) 

S-+00 


lim  v(t) 


(5) 


Thus  the  final  value  of  the  integrated  output  is  independent  of  the 
characteristics  of  the  (low-pass)  filter  and  no  account  need  be  taken  of  it  when 
calibrating  by  the  velocity  method  other  than  to  ensure  that  the  signal  has 
settled  to  its  final  value  before  shutting  off  the  integration. 

By  way  of  example.  Fig  2 shows  the  response  of  an  accelerometer  with  filter 
of  f = 28  kHz  and  £ » 0.01  to  a $ sin  pulse  of  unit  acceleration  and  of 
length  100  us  . The  sensitivity  of  the  accelerometer  is  assumed  to  be  unity, 
and  for  these  conditions  it  is  seen  that  if  the  peak  output  were  recorded  on  a 
peak-reading  meter  the  value  recorded  would  have  an  error  of  about  14%.  Moreover, 
a considerable  period  would  have  to  elapse  before  the  velocity  output  attained  a • 
steady  value. 

Fig  3 shows  the  error  in  measuring  peak  acceleration  for  various  values  of 

u)  (in  terms  of  T ) and  t , whilst  Fig  4a  gives  the  approximate  times  for  the 
n 

integrated  output  to  settle  within  1%  of  its  final  value.  It  is  evident  from 
Fig  3 that  the  error  is  very  small  for  values  of  (io^TJ/tt  >4  if  C = 1//2 
(Butterworth  filter). 

3 RESPONSE  TO  SIN2  PULSE 

Whereas  the  ^ sin  pulse  is  a useful  shock  waveform  for  calibration 

purposes,  under  severe  impacts  the  sharp  corners  at  the  beginning  and  end  of  the 

waveform  can  cause  excessive  excitation  of  the  transducer  resonance,  or  other 

resonances  associated  with  the  transducer  mounting.  Thus,  in  calibration,  it  is 

usual  to  place  a wad  of  felt  or  other  crushable  material  on  one  of  the  impacting 

2 

surfaces.  This  gives  a shock  waveform  not  unlike  a sin  pulse  (Fig  lb). 

2 

For  a sin  pulse  the  acceleration  is 


a(t)  = a sin  wt 

0 < t 4 tt/co 

(6) 

- 0 

t > tt/u)  . 

(6a) 

From  Appendix  B,  the  output  of  the  accelerometer  and  filter  is 


8 


y(t)  = g(t)  - g(t  - ^)H(t  - J) 


Zka^A  + B cos  2(ot  + sin  2<ot 
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If  the  accelerometer  and  filter  are  followed  by  an  integrator,  the  velocity 


output  is 


v(t)  = g(t)  - g(t  - >(t  - 


2ka  At  + B + C cos  2oit  + •$-  sin  2oot 
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The  final  value  of  the  velocity  output  is  again  independent  of  the  filter 
parameters  and  is  given  by 


lim  v (t)  = a . (9) 

t-wo 

Time  histories  of  the  acceleration  and  velocity  outputs  for  T = 100  ps  , 
f ■ 28  kHz  and  c,  = 0.01  are  shown  in  Fig  5.  Comparison  with  the  response  to 
a J sin  pulse  (Fig  2)  shows  that  both  the  ringing  and  peak  error  are  consider- 
ably smaller. 

The  error  in  measuring  peak  acceleration,  and  the  time  taken  for  the  inte- 
grated output  to  settle  to  within  1%  of  its  final  value,  are  shown  in  Figs  6 and 

4b  respectively.  The  error  is  small  for  values  of  (cu  T)/tt  > 4 if  c,  = 0.67  . 

n 

This  value  of  damping  ratio  is  close  to  that  of  the  so-called  ’optimum  damping’ 
value  used  in  transducers. 
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4 CALCULATION  OF  ACCELEROMETER  NON-LINEARITY  WHEN  CALIBRATING  BY  THE  VELOCITY 
METHOD 

Piezoelectric  sensors  in  general  display  a high  degree  of  linearity  and  so 
have  a very  wide  dynamic  range.  Nevertheless,  when  called  upon  to  measure  very 
large  accelerations  non-linearity  may  become  significant  in  transducers  with 
ceraaic  elements. 

In  practice,  the  sensitivity  tends  to  increase  linearly  with  acceleration 
input  so  that  it  has  the  form 

k'  = k(l  + Aa(t))  (10) 

where  k'  = true  sensitivity 

k = low-level  sensitivity 

and  A = fractional  change  of  sensitivity  per  unit  input  acceleration. 

Fig  7 shows  a sketch  of  the  sensitivity  and  input-output  relationship  with 
the  non-linearity  enhanced  for  clarity.  Clearly,  non-linearity  can  cause  pro- 
blems when  calibrating  by  the  velocity  method  since  point-by-point  calibration 
cannot  be  performed.  However,  if  A is  constent  (often  a fair  assumption)  then 
its  value  can  be  calculated  if  the  shape  and  duration  of  the  pulse  is  known. 

4. I Calculation  of  A when  calibrating  with  j sin  pulse 

For  a \ sin  pulse,  it  has  already  been  shown  (equation  (5))  that  the  final 
value  of  the  integrated  output  is,  with  integrator  of  transfer  function 
7 (gain  = |i|)  , 


v(t)  = 2 — a . 

t~>«>  11 


If,  however,  the  gain  is  k'  , the  final  output  will  be,  from  equations  (1) 
and  (10), 


v(t)  = a I k'  sin  cut  dt 

t-Hx.  I 


ka  / (sin  wt  + Aa  sin  wt)dt 
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The  fractional  error,  caused  by  non-linearity,  in  the  final  value  of  the 
integrated  output  is  therefore,  from  equations  (II)  and  (12), 


(13) 


where  E = fractional  error  in  velocity  measurement  caused  by  non-linearity. 

To  gain  an  estimate  of  the  value  of  X from  a velocity  calibration,  the 
output  recorded  at  one  or  more  high  acceleration  levels  must  be  related  to  the 
low-level  sensitivity  k obtained,  for  example,  by  calibrating  on  a vibrator. 
Equating  the  final  value  of  the  integrated  output  with  the  measured  velocity 
change  gives,  from  equations  (11)  and  (12) 


with 


v ( t ) = kTa  (—  + 

t-».  11  2 


Avtt 

~2T"" 


where 


Av  = measured  velocity  change. 
Hence 


X 


8t/v (t)  - 

\t-*=° 


kAv' 


2 2 
7i  k(Av) 


(14) 


. . 2 
4.2  Calculation  of  X when  calibrating  with  sin  pulse 

. 2 

For  a sin  pulse,  the  final  value  of  the  integrated  output  is,  from 
equation  (9), 


v(t) 

t-+« 


kTa 

2 


(13) 


If  the  gain  is  k’  , the  final  output  is,  from  equations  (6)  and  (10) 


12 


tt/oj 

.2  « . 4 

v(t)  = ka  j (sin  wt  + Aa  sin  wt)dt 

(;•+<» 


/ 


C 3Xa-> 

l1  + -H  • 


(16) 


The  fractional  error  is  therefore,  from  equations  (15)  and  (16),  given  by 


E = 


3Xa 


(17) 


From  equations  (15)  and  (16), 


v(t) 


- 


with 


Hence 


2Av 

T 


A = 


2l/v(t)  - kAv^ 


3k (Av) ‘ 


08) 


Calculation  of  A in  either  case  has  involved  a knowledge  of  the  pulse 

shape  and  period.  Whilst  the  latter  may  be  measured  with  fair  accuracy,  the 

former  must  be  a matter  of  qualitative  assessment.  Comparing  equations  (14)  and 

2 

(18)  shows  that  the  values  of  A for  i sin  and  sin  pulses  have  coefficients 
2 

of  8/ (tt  ) and  2/3  respectively.  Hence  the  value  of  A could  be  in  error  by 
as  much  as  20%  were  the  pulse  shape  wrongly  assessed.  However,  since  A is,  or 
should  be,  small,  errors  in  its  measurement  should  not  be  of  great  significance. 
Typically,  the  value  of  A is  such  that  k'  at  the  upper  limit  of  range  is  5% 
higher  than  k . 

The  above  treatment  assumes  that  A is  constant  and  independent  of 
acceleration.  The  validity  of  this  assumption  can,  of  course,  be  confirmed  by 
performing  velocity  calibrations  at  a number  of  acceleration  levels. 

5 EXAMPLES  RELATING  TO  THE  AIRGUN  USED  FOR  SHOCK  ACCELEROMETER  CALIBRATION 

To  illustrate  some  applications  of  the  data  derived  in  previous  sections, 

typical  examples  relating  to  the  calibration  of  accelerometers  in  RAE  with  an 

airgun  facility  are  considered  here.  The  airgun  produces  an  acceleration  pulse  of 
6 — 2 

up  to  10  ms  with  a period  of  about  I00  ps. 
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5. 1 Relative-motion  error  when  calibrating  by  the  comparison  method 

In  calibrating  by  the  comparison  method,  the  test  accelerometer  is  mounted 
on  the  end  face  of  a special  reference  accelerometer  which  has  been  previously 
calibrated.  The  peak  outputs  of  the  test  and  reference  accelerometers  are 
recorded  on  peak-reading  voltmeters  and  the  readings  compared.  This  method  is 
both  simple  and  reliable  provided  that  the  test  accelerometer  experiences  the 
same  input  acceleration  as  the  sensor  in  the  reference.  Errors  may  arise,  how- 
ever, because  there  is  in  practice  some  relative  motion  between  the  test 
accelerometer  and  sensor. 

It  has  been  deduced  from  the  manufacturer's  data,  and  confirmed  by 
frequency-response  measurements  on  a vibrator,  that  the  mechanical  couplings 
between  the  accelerometer  and  reference  transducer  normally  used  may  be  represen- 
ted by  a second-order  system  having  a natural  frequency  of  29  kHz  and  damping 
ratio  of  0.49.  Using  these  values  in  equations  (3)-(3e)  and  (7)-(7f)  shows  that 

the  peak  acceleration  applied  to  the  test  accelerometer  is  1.4%  higher  than  that 

.2 

experienced  by  the  reference  for  a ^ sin  pulse  and  3.0%  higher  for  a sin 
pulse.  Thus,  although  a sin  pulse  produces  less  ringing  than  a i sin  pulse, 
the  relative-motion  error  may  be  greater.  Nevertheless,  relative-motion  errors 
are,  in  general,  small. 

5. 2 Error  caused  by  accelerometer  resonance 

Shock  accelerometers  typically  have  an  undamped  natural  frequency  of  about 

100  kHz  with  small  damping.  Thus,  ringing  will  occur  and  the  amplitude  of  the 

peak  registered  will  be  in  error.  Fig  8 shows  the  error  in  measuring  the  peak 

2 

acceleration  of  | sin  and  sin  pulses  of  100  us  duration  for  transducer 

natural  frequencies  ranging  from  50-150  kHz,  assuming  a damping  ratio  of  0.01. 

Since  the  damping  ratio  is  very  small,  the  greater  ringing  excited  by  the  i sin 

2 

pulse  produces  larger  errors  than  those  with  a sin  pulse.  The  effects  of 
accelerometer  resonance  can  often  be  reduced  by  careful  use  of  electrical  filters. 

6 CONCLUSIONS 

It  is  hoped  that  the  data  presented  in  this  Memorandum  may  afford  the 
reader  a better  picture  of  some  aspects  of  the  performance  of  shock  accelerometers 
and  provide  a means  of  quantifying  some  of  the  errors  which  can  arise  during 
calibration.  Moreover,  the  data  on  filtered  response  should  be  useful  in 
applications  where  mechanical  filters  are  deliberately  placed  in  the  measurement 
chain. 

More  specifically,  it  has  been  shown  that:- 


14 


(1)  The  best  filter  for  a j sin  pulse  is  one  with  a damping  ratio  of  \//l  , 

2 . 

whilst  for  sin  pulses  a damping  ratio  of  about  0.67  is  the  optimum. 

(2)  The  parameters  of  the  filter  do  not  affect  the  final  value  of  the  integral 
of  the  accelerometer  output. 

(3)  The  error  associated  with  relative  motion  when  calibrating  by  the 
comparison  method  is  unlikely  to  exceed  3%  with  the  accelerometer  and 
reference  transducer  presently  used  in  RAE. 

(4)  It  is  fairly  easy,  in  principle,  to  measure  an  approximate  value  of 
accelerometer  non-linearity  when  calibrating  by  the  velocity  method, 
provided  that  the  accelerometer  sensitivity  increases  or  decreases 
linearly  with  acceleration. 

Clearly,  caution  is  necessary  in  dealing  with  practical  situations  where 
mechanical  filters  may  display  characteristics  other  than  those  of  a second-order 
system,  and  where  the  shock  differs  significantly  from  the  two  ideal  pulse  shapes 
considered. 
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Appendix  A 

RESPONSE  TO  j SIN  PULSE 

For  an  input  of  a(t)  = a sin  <ot  , 0 < t { - 


= 0 , 


IT 

t > — 
(0 


its  Laplace  transform  is 


= ajsin  ut  + j^sin  <*»(«:-  J)j  Jh  (t  - | 


a(s)  = a 


w__  T + e~ (sir) /o)"l 

+ ID  *" 


2 

S + ID 


If  a second-order  filter  of  transfer  function 


10 


h(s) 


n 


2 2 
s + 2cto  s + to 
n n 


is  placed  before,  or  after,  the  accelerometer  of  sensitivity  k , the  transform 
of  the  output  is 


y(s)  = ka 


ka 


r + -(sto/io] 

(S2  + (02)  (s2  + 2C(0  s + (02)  L J 

K n nJ 


(A—  1 ) 


As  + B Cs  + D 


2 2 2 2 
s + (o  s + 2 £<o  s + w 
n n 


[l  ♦ e"(sit)/o)] 


= ka 


As  ->  B , C<8  + * ° ' CWn 

s2  + «i2  (s  » cu  )2  ♦ u2(l  - C2) 
n n 


[>  + e-(sTt)/u.J 


with  the  coefficients  of  the  partial  fractions  being 
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0)  0) 

“ 
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- 

nJ 

[r  2 2-j  2 2 2 2] 

i Mto  - to^J  + 4c  to  to  I 


Inverting  the  transform  gives 

a f E . B « 


-£to  t 


a(t)  *=  ka  j^A  cos  tot  + — sin  <ot  + Ce  n cos  to  (1  - c2)^t 

D - CCto  -Ct o t . 

+ jT  e n sin  to  (1  - c/ t + g(t--)H(t--) 

to  (1  - C V w u 

n 

L — 


(A- 2) 


with  g(t  - Jj  being  the  first  term  in  the  square  brackets  with  (t  - — ) 
substituted  for  t , and  H(t  - being  Heaviside’s  shifting  function. 

If  an  integrator  is  added  to  the  system,  equation  (A-l)  becomes 


u j- 

y(s)  = ka  j - -—2 J 

s(s  + to  ) fs  + 2cto  s + to2)  ^ 
n n' 

This  may  be  rewritten  as 


I + e 


■(sn)/io  J 


(A- 3) 


n n' 


y(s)  - ka 


A Bs  + C D(s  + C<*>)  + E — Dcto  «- 

* ? * -T-H  * ~§ 5 f-  fl  * e-(s,)/"l 

s + to  (s  + ^ton)  + u/j(l  - c ) L J 


with  the  coefficients  being 
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Inverting  the  transform  gives  the  velocity  output  of 


v(t)  * ka 


C “C<*>  C 9 i 

A + B cos  (Dt  + — sin  cot  + De  n cos  w (I  - t 
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Appendix  B 

RESPONSE  TO  SIN2  PULSE 


. 2 


For  a sin  pulse,  the  input  is 


a(t) 


" • 2 „ 
a sxn  wt  , 


0 < t v<  - 
u> 

t > — 

(0 


Since 


a sin  uit  = 


a ,,  2 . 

(1  ~ cos  wt)  , 


a(t) 


Hence 


f j[l  - COS  2wt]  - [l  - cos  2„(t  - 1)][h(c  - 1)]| 
!<8)  ■ 

- 3 4a)2  f,  _-(S7T)/w"l 

2 2 2 e • 
s(s^  + 4w  ) L J 


With  the  addition  of  a second-order  filter  the  transform  of  the  output 
becomes 


y(s) 


2a 


2 2 
w w 

n 


,2  . 2.  , 2 „ 2'. 
s(s  + 4w  ) ( s + 2i;w  s + w 
' n nJ 


e 


(S7r)/w 


Splitting  this  into  partial  fractions,  as  performed  in  Appendix  A,  and 
inverting  the  transform  leads  to  the  output-time  relationship  described  by 
equation  (7).  A similar  process  leads  to  equation  (8)  describing  the  velocity 
relationship. 
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